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ABSTRACT 

We present an analysis of X-ray colour maps of the cores of clusters of galaxies, formed from 
the ratios of counts in different X-ray bands. Our technique groups pixels lying between con- 
tours in an adaptively-smoothed image of a cluster. We select the contour levels to minimize 
the uncertainties in the colour ratios, whilst preserving the structure of the object. We extend 
the work of Allen & Fabian (1997) by investigating the spatial distributions of cooling gas 
and absorbing material in cluster cores. Their sample is almost doubled: we analyse archive 
ROSAT PSPC data for 33 clusters from the sample of the 55 brightest X-ray clusters in the 
sky. Many of our clusters contain strong cooling flows. We present colour maps of a sample of 
the clusters, in addition to adaptively-smoothed images in different bands. Most of the cooling 
flow clusters display little substructure, unlike several of the non-cooling-flow clusters. 

We fitted an isothermal plasma model with galactic absorption and constant metallicity 
to the mid-over-high energy colours in our clusters. Those clusters with known strong cooling 
flows have inner contours which fit a significantly lower temperature than the outer contours. 
Clusters in the sample without strong cooling flows show no significant temperature variation. 
The inclusion of a metallicity gradient alone was not sufficient to explain the observations. A 
cooling flow component plus a constant temperature phase did account for the colour profiles 
in clusters with known strong cooling flow components. We also had to increase the levels of 
absorbing material to fit the low-over-high colours at the cluster centres. Our results provide 
more evidence that cooling flows accumulate absorbing material. No evidence for increased 
absorption was found for the non-cooling-flow clusters. 

Key words: galaxies: clusters: general - cooling flows - intergalactic medium - X-rays: 
galaxies. 



1 INTRODUCTION 

Clusters of galaxies are among the largest dynamical objects in the 
universe. They contain large quantities of virilised gas (the intr- 
acluster medium, or ICM) at temperatures of the order of 10 7 - 
10 8 K, heated by their initial gravitational collapse. Due to this 
gas they are powerful emitters of X-rays from bremsstrahlung and 
line emission. The gas masses of rich clusters exceed 10 14 M@, and 
they can emit over 10 45 erg s _1 in X-rays. 

The radiative cooling time, i CO oi, is shortest in the densest re- 
gions, at the centre of a cluster, and is often found to be less than 
the plausible age of the cluster, i.e. the time since its last major 
interaction, which is a significant fraction of the Hubble time. As 
the gas cools, its temperature decreases, but its density must rise to 
maintain hydrostatic equilibrium. This is achieved by the gas above 
flowing in. This is known as a cooling flow (See Fabian 1994 for a 
review). 

Observations have shown that cooling flows are inhomoge- 
neous. The mass drops out of the flow into a form which does not 
emit X-rays. The exact destination and final form of this gas has 
not been determined, but magnetic fields are probably important in 
determining, and constraining, the structure of cooled gas. 



The rate of the mass of gas cooling out at a certain radius is 
known as the mass deposition rate, M(r). This value can total from 
10-1000 M Q yr _1 . The mass deposition rate can be estimated by 
assuming that all the luminosity within the cooling radius, where 
the cooling time is less than the age of the cluster, is due to the 
cooling gas radiation, plus the PdV work done on the gas as it en- 
ters that radius. The relationship between the mass deposition rate 
and radius is found to be roughly linear, as the surface brightness 
profiles are not as centrally peaked as would be expected for a ho- 
mogeneous cooling flow. 

Allen & Fabian (1997) undertook a study of the spatial distri- 
butions of cooling gas and intrinsic X-ray absorbing material in a 
sample of nearby clusters, using X-ray colour profiles, formed from 
the ratio of flux in selected bands, using data from the ROSAT Po- 
sition Sensitive Proportional Counter (PSPC). They found that the 
profiles indicated that there were significant central concentrations 
of cooling gas in their cooling flow clusters, becoming approxi- 
mately isothermal at larger radii. Their profiles also showed large 
levels of X-ray absorption in the cluster cores, increasing with de- 
creasing radius. 

The bands they examined are listed in Table [ll which they 
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Table 1. Definition of the energy bands (Allen & Fabian 1997) 
Band Energy range (keV) 



Ax 
M 

B 
C 
D 
F 



0.10-0.20 
0.20-0.40 
0.41-0.79 
0.80-1.39 
1.40-2.00 
0.41-2.00 



combined to form two X-ray colour ratios, B/D and C/D. To 
analyse the observations, they calculated theoretical curves for the 
colours as functions of temperature and metallicity, with a uniform 
screen to model the absorbing column density. They modelled both 
colours with isothermal and cooling gas models. In Section ^ we 
reevaluate these functions using more recent plasma models. 

To demonstrate the presence of distributed cooling gas in the 
central regions of cooling flow clusters, they used the C/D pro- 
file, and the B / D ratio to demonstrate the presence of intrinsic X- 
ray absorbing gas. They found the absorbing material only partially 
covers the X-ray emitting region. 

In this paper, we extend the analysis of Allen & Fabian (1997). 
Instead of making a simple radial profile of each cluster, we have 
developed a technique to take account of the two-dimensional 
structure of the image. Simple binning works in areas of high flux, 
where the errors created in forming ratios of counts are small. In 
low flux regions, the number of pixels required to bin destroys most 
of the spatial information in the image. Our method attempts to 
make statistically useful X-ray colours, whilst preserving spatial 
resolution. 



2 OBSERVATIONS AND COLOUR MAP ANALYSIS 

We analysed archival ROSAT PSPC data for a sample of 33 clus- 
ters. The clusters we examined are shown in Table 0. The sample 
was mainly chosen to include known strong cooling flow clusters 
(from Peres et al. 1998), with a few non-cooling-flow clusters for 
comparison. 

We used the bands of Allen & Fabian (1997) in our analysis. 
These bands are optimised for the study of cooling gas. This pre- 
cluded the use of the extended object analysis package of Snowden 
& Kuntz (1998), although we made use of some modified versions 
of a few of their routines. 

The full analysis procedure for each cluster is as follows. 

(i) We removed time periods from the observation where the 
Master Veto count rate averaged over ~ 30 s was greater than 220 
counts s _1 (Snowden et al. 1994). 

(ii) A modified version of Snowden's CAST_DATA routine was 
used to create 512 x 512 pixel images of the cluster, at a pixel size 
of approximately 15 x 15 arcsec 2 . We created images of the cluster 
in the B, C and D bands from Table 

(iii) We corrected the B, C and D band images for background 
counts (Section 2.2). 



(iv) An adoptively- smoothed image (Section 2. 1 ) was made of 



the central region of the cluster, in the high energy D band. We 
created an image of size 64 x 64, or occasionally 128 x 128 pixels, 
depending on the extent of the cluster. 

(v) We created linear contours in the adaptively-smoothed im- 
age, by simply taking the difference between the maximum and 
minimum count rates, and dividing by the number of contours. 



(vi) The counts for those pixels between each two neighbouring 
contours was added, in each of the background corrected B, C and 
D images. 

(vii) We used the count totals to generate average B/D and 
C / D colours for each contour. The error on each result is: 



2_/pix 



(TV* + C?) 



pixels C^i 

E • i [nY 

•i-— 'pixels V * 



+ C} 



^^-^ pixels C^i 



C} 



(1) 



where r is the X/Y average ratio, and o r is the error on the ra- 
tio, is the uncorrected number of photons in pixel i in band X, 
and Cj X is the background correction for pixel i in band X. This ex- 
pression does not take account of systematic background uncertain- 
ties. We ignored pixels that were suspected to contain point sources 
(Section 2.2) outside the core. 



(viii) A Monte Carlo procedure was used to change the contour 
levels in the adaptively-smoothed image to minimize the total error 
squared in the contours for both colours, or the value 



(B/D), J 



+ * 



(C/ 



D)J 2 }- 



(2) 



The result of moving the contour levels to minimize x creates con- 
tour levels which contain similar numbers of counts and signal to 
noise ratios. The procedure produces colour ratios with roughly 
equal fractional uncertainties for each contour. 

We used six contours to process the cluster images. This was 
a good compromise between showing detail and keeping statistical 
errors low. The results were stable to small changes in the number 
of contours. 



2.1 Adaptive smoothing 

Adaptive smoothing is a technique to smooth an image and re- 
tain significant small-scale structure (Ebeling, White & Rangarajan 
2000). Smoothing is required because it is often the case, especially 
in X-ray data, that there are few photons per pixel. However, con- 
volving an image with a Gaussian kernel leads to broadening of 
features and loss of small-scale structure. An algorithm is required 
which retains information on all scales. 

An adaptive kernel smoothing (AKS) algorithm applies a ker- 
nel with a varying size to an image. A part of the image with few 
counts is smoothed with a large kernel, whilst a part with many 
counts is smoothed with a small kernel. The algorithm decides 
what is the 'natural' smoothing scale for a part of an image. AKS 
algorithms differ on how they choose an appropriate kernel size. 
We used the ASMOOTH algorithm of Ebeling, White & Rangara- 
jan (2000) to adaptively smooth our cluster images. ASMOOTH is 
unusual in that it uses the local signal-to-noise ratio to adjust the 
size of the kernel. The minimum signal-to-noise ratio of the part of 
the image smoothed by the kernel is a parameter to the algorithm. 
A level of 6-a seemed to produce the best results; it appeared not 
to generate false detail from spurious pixels. We used a Gaussian 
rather than top-hat adaptive kernel. 

Other methods for revealing structure in data include wavelet 
analysis (Starck & Pierre 1998). Different methods have other ad- 
vantages, but the ASMOOTH algorithm is particularly simple and 
quick to perform. 
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2.2 Background removal 

The raw data contains counts not from the observed source. These 
background counts, sources of which are listed below, must be re- 
moved. 

(i) The X-ray background, which comes through the aperture of 
the telescope, and is observed with the same vignetting response 
of the telescope as the object. This includes solar and extragalactic 
background emission. 

(ii) The particle background (Plucinsky et al. 1993) is diffuse 
and is caused by high energy cosmic rays. The flux of the particle 
background is approximately constant across the detector. 

These components vary in strength with X-ray energy, so they must 
be removed to make X-ray colours. 

To estimate the background we took the average count from 
the 30-35 arcmin region, beyond the 'ribs' in the field of the PSPC. 
From this count we subtracted an estimate of the particle back- 
ground from Plucinsky et al. (1993), of 4 x 1CP 6 counts pixel -1 
s -1 arcmin -2 keV -1 . We corrected the result for vignetting. The 
corrected count and the particle background was removed from all 
pixels. 

Accounting for the total background reduces the scatter of the 
observed data points. It also hardens the emission from the outer- 
most parts of clusters which emit few photons. 

Three clusters, A2142, A1795 and A2029, contained obvious 
point sources in their fields outside the cluster core. We removed 
the point sources using a mask which excluded the affected region 
around them. Those pixels inside the mask were not included in the 
sum of counts inside a contour. 



X-ray colour maps of the cores of galaxy clusters 
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Table 2. The sample of clusters examined, z is the redshift taken from the NASA/IPAC Extragalactic Database (NED), and 
the values of M, the mass deposition rate, are PSPC observation derived values taken from Peres et al. (1998), as are the 
central cooling time values, t coo i. The sequence identification and exposure time give PSPC observation details. S(B/D) 
and S(C/D) are the significance of a colour gradient, larger values imply a larger colour gradient (Section [?]). The clusters 
are listed in cooling flow flux order, M / z 2 . 



Cluster 


Sequence ID 


Exposure/s 


z 


M/(M Q yr" 1 ) 


tcool/Gyr 


S{B/D) 


S(C/D) 


Virgo 


rp800187n00 


10536 


0037 


oy_ 9 


2+ 00 




21 


Perseus 


rp800186n00 


4711 


U.UIOJ 


c^fi+33 
OOu 24 


n q+ 00 


S Q 
j, j 


Q A 
0.4- 


2A 0335+096 


rp800083nOO 


10220 


n 0^4Q 


qqc:+32 
OZJ -43 


n q+ 00 




7 1 


Centaurus 


rp800192n00 


7793 


0109 


on + 10 

5 


8+ 00 


1 9 


14 


A2199 


rp800644n00 


40QQQ 


030 


8 


1 q+00 
i- 9 -o.i 


7 8 


1 


Ophiuchus 


rp800279n00 


3932 


028 


1 27+ 48 


o n +0.4 


-1 2 


n t 

U. J 


Klemola 44 


rp800354n00 




n 098^ 


8 '-19 


9 + O.6 

z - o -0.2 


-1 6 


n 1 


A2052 


rp80O275n0O 


6211 


0348 


i nc+26 


2 c ;+ - 1 


2 1 


d (S 
t.D 


A262 


rp800254n00 


GAGA 
OOOO 


0164 


27+ 4 
z '-3 


1 c i+ - 1 


2 2 


J. 1 


PKS 0745-191 rp800623n00 


1 047^ 


n 1078 


1 O^S+ 116 
!Ud8_ 6g 


9 9+O.I 

z - z -o.i 


z.o 


1 7 
Z. / 


A 1060 


rp800200n00 


15764 


0124 


1C+5 


4 7+0- 4 
4 -'-0.3 


3 


(1 1 

U. I 


A1795 


rp800105n00 


jDi / J 


0627 


381 + 41 


l- y -0.1 




S Q 


Hydra-A 


rp800318n00 


1 0J70 


O OS99 


9fiA+^l 
zb4 -60 


9 n +o.o 


■3,0 


c c 
J.J 


A2029 


rp800249n00 


12542 


0767 


73 


2 0+0.1 


4 4 


S 1 

j .z 


A496 


rp800024n00 


GGS7 


033 


qc + 13 

yo_ 12 


1 8+ - 1 


1 4 


J.o 


MKW3 


rp800128n00 


9984 


O 044Q 


17K+14 

1 '°-46 


o n +o.i 


n 1 
-u. 1 


n r 

U.O 


A478 


rp800193n00 


21969 


O OSS? 


fi1 fi+ 63 
oio_ 76 


2 8+ - 1 


-U.o 


2 1 


A85 


rp800250n00 


10238 


0521 


lyo -52 


2 4+ - 1 


{ 5 


2.9 


A3112 


rp800302n00 




0746 


Q7«+80 
■J'O-61 


1 Q+0-1 

i- y -o.i 


1 4 


J.O 


Cygnus-A 


rp80()622n00 


9442 


057 


244+ 26 


2 fi+ 01 


U.D 


1 (S 

1 ,0 


AWM7 


rp800168n00 


13088 


0172 


41+6 


1 q+0- 2 


-0.4 


2 2 


A4059 


rp800175n00 


5439 


0.0478 




„ 4+O.6 


0.9 


1.9 


A3571 


rp800287n00 


6062 


0.0391 


79+56 
' z -31 


°- s -i.o 


6.7 


3.4 


A2597 


rp800112n00 


7163 


0.0824 


271±« 


9 o+O.l 

z -°-o.i 


-0.3 


1.4 


A644 


rp800379n00 


10246 


0.0704 


189±^ 6 


6.8ir 4 


2.8 


1.4 


A2204 


rp800281n00 


5357 


0.1523 


ocn + 127 


3.i±g:i 


2.8 


3.4 


A2142 


rpl50084n00 


7734 


0.0899 


OJU_ 133 


c 9+0.4 


1.0 


1.0 


A3558 


rp800076n00 


29490 


0.0478 


40+f 


io.2l°; 3 2 


-0.2 


1.2 


A401 


rp800235n0O 


7457 


0.0743 




10.6+*° 


-0.6 


-1.3 


Coma 


rp80()005n00 


21140 


0.0232 


n +1 
u -o 


i7.7il;l 


1.7 


0.9 


A754 


rp800550n00 


8156 


0.0542 


n+29 
U -0 


15.0+2 2 


1.6 


2.5 


A2256 


rp800162a01 


4747 


0.0581 


n+14 
U -0 


15-Ojlg 


0.2 


0.9 


A119 


rp800251n00 


15197 


0.044 


n +2 


1 q 9+12-2 


-0.9 


-1.8 
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Table 3. Mean distance of the pixels in each contour from the cluster centre. 
Distances are measured in pixels (15 arcsec). The contours are numbered 
from 1 (inner) to 6 (outer). Clusters marked with an asterix (*) were ex- 
amined using a 128 X 128 pixel adaptively smoothed image, rather than a 
64 X 64 image. 



Cluster 


1 


2 


3 


4 


5 


6 


2A ujjj+uyo 


2. 1 


1 c 
3.5 


j.y 


5. / 


IT Q 

13.0 


9C "7 

25. / 


A 1 AAA* 

AlUOU^ 


A o 

o.y 


13.2 


im 

iy.z 


nA ^ 
2o. J 


17 A 

3 /.U 


AA 7 
OU. / 


Al iy^ 


"7 c 

/.J 


AAA 
14.0 


20.5 


11 .J 


n a 
3/.U 


Af\ A 

00.4 


A 1 "70^ 

Al tyj 


9 A 

2.0 


4.2 


A A 
0.0 


Q A 

y.o 


1/1 T 
14. / 


25. y 


a imn 
A2U2y 


2. / 


A 1 

4.3 


A A 

0.4 


y.j 


13.7 


25. J 


A2UJ2 


3.U 


j.U 


7 O 

/.y 


Ifl c 

1U. J 


1 A Q 

14. y 


70 

25. y 


A2142 


3.1 


5.2 


8.3 


12.2 


18.3 


29.5 


ATI QO 


1 A 
3.0 


7 A 


Ifl c 

1U.J 


1/1/1 
14.4 


OA 1 

2U.3 


11 1 

31.1 


A22U4 


2.0 


2.y 


1 A 

3.0 


c c 

J. J 


5.0 


ia n 

20. y 


AZZJO 


"7 /I 
1 A 


11 Q 

i i.y 


1 A 7 

lo. / 


OA Q 
2U.5 


9< A 
2J.D 


<A 1 

DO. 1 


A 

A2jy / 


O /I 

2.4 


9 o 

2.y 


1 7 
3. / 


J. 3 


C 1 

5. 1 


97 A 

2/.U 


AZoZ 


A 1 
0.3 


11 o 
11.5 


lO 1 

15. 1 


9/1 Q 
24.5 


11 A 

33.4 


<a 1 

jy. 1 


Aim 
A31 12 


3.3 


'J 7 
3. / 


^ 9 
J. 2 


Q ft 

5.U 


111 
13.1 


70 c 
25. J 


Ajjjo p 


7 1 
/.I 


1 9 ^ 
12. J 


1 IX) 


9 1 A 
21 .U 


97 A 
Z/.u 


S,A 8 
JO. 5 


A 1^7 1 * 

A3 J / 1 * 


c c 

J.J 


1 A 9 
1U.2 


AAA 
14.4 


21. / 


19 1 

32. 1 


J5.U 


A /1fl1 

A4U1 


1 A 
3.U 


0.4 


y. / 


12.5 


1 A o 

lo.y 


90 A 

2y.o 


A /I A^O 


i i 
3.3 


^ 7 
J. 1 


O.J 


11 c 
1 l.J 


1 /.j 


mi 

3U.3 


rV^t / 


9 4 
Z.4- 


3.D 


J. 1 


/ .D 


1 1 9 
1 l.Z 


97 8 
Z / .0 


A /I OA 

A4yo 


2.5 


J. / 


O A 

y.4 


13. y 


9A 1 
2U. 1 


1 1 A 

31. U 


A AA /I 
A044 


3.5 


A 1 
0. 1 


Q 7 
5. / 


11 o 

11.5 


1 A A 
10.4 


TO A 

zy.u 


A754* 


12.0 


25.7 


27.7 


30.7 


36.6 


57.2 


A85 


2.9 


5.9 


10.0 


14.0 


19.5 


30.8 


AWM7* 


7.9 


14.4 


21.1 


27.7 


36.5 


60.0 


Centaurus* 


5.7 


12.3 


18.5 


27.2 


39.2 


60.9 


Coma* 


14.0 


22.5 


29.4 


36.7 


46.6 


63.0 


Cygnus A 


4.6 


5.3 


6.2 


8.8 


16.5 


29.5 


Hydra A 


2.6 


3.8 


6.7 


9.6 


13.4 


28.4 


KLEM44 


4.3 


6.0 


8.9 


12.4 


17.8 


29.8 


MKW3 


2.4 


4.7 


7.0 


9.2 


12.9 


28.2 


Ophiuchus 


7.3 


12.2 


15.2 


19.1 


24.4 


31.3 


Perseus* 


4.2 


8.4 


13.5 


22.2 


34.8 


59.3 


PKS 0745-191 


2.4 


3.1 


4.0 


6.2 


9.6 


27.3 


Virgo* 


7.4 


14.2 


21.7 


30.2 


41.1 


62.2 



© 0000 RAS, MNRAS 000, 000-000 



6 J.S. Sanders, A. C. Fabian and S. W. Allen 



3 RESULTS 

We created adaptively-smoothed images of the sample of clusters in 
the B, C and D bands, using an ASMOOTH minimum significance 
of 6-ct. Fig. |l| ^ and ^ show the images for the Perseus, Virgo and 
A2199 clusters, respectively. ROSAT HRI images of the centres of 
the Perseus cluster, NGC 1275, and the Virgo cluster, M87, have 
been examined by Bohringer et al. (1993, 1995), respectively. They 
found the thermal plasma was displaced by the radio lobes of NGC 
1275 at the centre of the Perseus cluster. 

Fig. g and | are colour contour maps of the Perseus and Virgo 
clusters. The maps show the average X-ray colours between each 
of six contours. The scale below each map shows the value and sta- 
tistical uncertainty of the colour of each contour. The extreme error 
bars on the softest and hardest points are not shown on the scales. 
Note that darker shades in these plots indicate softer emission. Ta- 
ble U shows the mean physical distance from the cluster centre of 
each contour in each cluster in pixels. It also identifies which clus- 
ters were examined using a 64 x 64 pixel image or a 128 x 128 
image. 

Fig. |^ shows a sample of six colour contour maps. Four of the 
clusters are non-cooling-flow clusters, showing interesting features. 
Two are cooling flow clusters, which make up most of our sample, 
and they reveal little substructure. The point source in A2142 is 
shown, but is not included in the contour statistics. 

Fig. ^ and || show the average B / D and C/D colours for the 
contours of each of the sample of clusters. The clusters are listed 
in decreasing cooling flow flux order (M / z 2 , where M is the mass 
deposition rate, and z is the redshift of the cluster) using the PSPC 
M values obtained by Peres et al. (1998). However, some of the 
clusters have highly uncertain mass deposition rates, so the order- 
ing is not definite. 

Table ^| shows a significance measure of the change in colour 
from the outside to the inner region for each cluster. For a cluster 
with contours of colour Xi, it is calculated using 



S(X) 



X e )/2] : 



(3) 



Xi is the innermost colour value, Xs and Xe are the two outermost 
colours, ax is the statistical uncertainty in each colour. The uncer- 
tainties are almost the same for each contour, due to our method 
of choosing contour levels (Section [2]). \/3/2 is a numerical fac- 
tor to make S(X) a measure of the number of 1-ct errors. We take 
the mean outer two colours to decrease the chance of background 
contamination. The background is a large fraction of the total num- 
ber of counts in the outer contour for several clusters, for example 
A3571, A2204, A3558 and A262. 

Fig. ^| shows the variation of the significance of the C/D 
colour, presumably due to the presence of cool gas, with the cool- 
ing flow flux, M/z 2 , for the sample of clusters. Fig. shows the 
significance plotted against central cooling time, i coo i. There is a 
clear anti-correlation between t coo i and the significance of a C/D 
colour gradient. Clusters with t coo i < 2 Gyr show colour gradients 
with S(C/D) > 2. t coo i is a measure of the quality of data for a 
cooling flow. Cooling flow clusters with short cooling times have 
enough counts to make their central bins small in area. 

Some of the results obtained are dependent on our choice of 
significance measure. If the cluster has a large soft core, as A3571, 
Hydra-A and A2204 seem to have in B / D, then using a measure 
based on X3 + X4 instead of X5 + Xe in equation (J3[> gives a much 
lower significance value. A significance based on X3 + X4 would 
be less affected by background noise, but throws away much of the 



o 

CV) 



o 2 



1000 



10* 10° 10° 

(M/z 2 )/(M yr-- 1 ) 



Figure 9. Plot of cooling flow flux, M / z 2 . against the significance measure 
of a colour gradient S(C/D) (Section H). The 2-a significance level is 
indicated by a dotted line. 
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Figure 10. Plot of the central cooling time, i coo i against the significance 
measure of a colour gradient S(C/D) (Section R|). The 2-a significance 
level is indicated by a dotted line. 



data. A more complex significance measure might be useful, but 
for those clusters where it might be important, noise is probably a 
more limiting factor. 
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Figure 4. Contour colour maps of the Perseus cluster. Areas between contours are coloured by their average X-ray colour. The points below the maps show 
the numerical colours, with their Poisson errors. The most central contours are the uppermost points. Darker areas indicate softer emission. 



Virgo in B/D Virgo in C/D 
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Figure 5. Contour colour maps of the Virgo cluster (see Fig. Q for a description). 
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Figure 7. Variation of the ratio B/D for the sample of clusters. Each point represents the colour of a contour. The uppermost points represent the innermost 
contours, numbered 1. 
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Figure 8. Variation of the ratio C/D for the sample of clusters. Each point represents the colour of a contour. The uppermost points represent the innermost 
contours, numbered 1. 
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Figure 11. Theoretical predictions for the B/D and C/D ratios for an isothermal gas at temperatures of 1-10 keV. The lines show steps of 1 keV in 
temperature. The lowest temperature at the highest ratio on the j/-axis. (a) B/D and (b) C/D as functions of metallicity at fixed column density (N H = 
10 20 atom cm -2 ), (c) B/D and (d) C/D as functions of column density at fixed metallicity (Z = 0.5Zq). 




Figure 12. Theoretical predictions for the B/D and C/D ratios for an gas cooling from temperatures of 1-10 keV to 0.001 keV. The lines show steps of 
1 keV in temperature. The lowest temperature at the highest ratio on the j/-axis. (a) B/D and (b) C/D as functions of metallicity at fixed column density 
(N H = 10 20 atom cm -2 ), (c) B/D and (d) C/D as functions of column density at fixed metallicity (Z = 0.5Zq). 



© 0000 RAS, MNRAS 000, 000-000 



X-ray colour maps of the cores of galaxy clusters 13 



(a) 



(b) 



(c) 



(d) 



Figure 13. Temperature profiles generated by fitting constant galactic obscuring column density and constant 0.3 Zq abundance to the C/D colour profile. 
The clusters are (a) A1060 (non-cooling-flow), (b) A2199, (c) Centaurus and (d) Virgo. The errors are 1-cr, generated from the errors in C/D. The contours 
are numbered from innermost to outermost. 
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Figure 14. Profiles showing how the fractional cooling flow component must increase to fit the C/D ratio assuming an isothermal temperature, constant 
galactic column density and constant 0.3 Zq abundance. The clusters are (a) A1060 (non-cooling-flow), (b) A2199, (c) Centaurus and (d) Virgo. The errors 
are 1-cr, generated from the errors in C/D. The contours are numbered from innermost to outermost. 
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Figure 15. Profiles showing how the intrinsic absorption must increase assuming an isothermal temperature, constant galactic column density and a constant 
0.3 Zq abundance. These were calculated by assuming the cooling flow components above, and calculating the extra required column density above galactic 
to explain the B/D ratio, assuming constant 0.3 Zq abundance. The errors are 1-cr, and were generated from the error in B/D, the error in the cooling flow 
component is ignored. The clusters are (a) A1060 (non-cooling-flow), (b) A2199, (c) Centaurus and (d) Virgo. The contours are numbered from innermost to 
outermost. 
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Table 4. Table showing the clusters and their calculated excess absorbtion, 
AN H (units 10 20 atom cm -2 ). Also shown is the used galactic column 
density, N H ga} (units 10 20 atomcm~ 2 ), and the significance of the excess 
absorption, S(AN H ). The clusters are listed in cooling flux order (M/z 2 ). 
The outer contour of A357 1 was ignored as it was suspected to have back- 
ground contamination. Ti so is the isothermal temperature calculated from 



the outer contours, which is limited to a maximum of 10 keV. 
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71 «„ / keV 
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Virgo 
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13.7 ±0.7 


19.0 


2.5 ± 0.1 
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14.9 


3.1 ± 1.1 
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2A 0335+096 
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22.7 ±2.4 
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4 COLOUR MODELS 

By comparing the B / D and C / D ratios with theoretical models, 
it is possible to estimate the effective temperature, metallicity and 
absorbing column density. We reevaluated the theoretical curves of 
Allen & Fabian (1997) using a more recent model of an isother- 
mal gas; the XSPEC 10 MEKAL model based on the calculations of 
Mewe and Kaastra with Fe L calculations by Liedahl (Mewe, Gro- 
nenschild & van den Oord 1985; Liedahl, Osterheld & Goldstein 
1995) with an absorbing PHABS screen. 

We iterated the calculation of the B / D and C / D colour ratios 
over the parameter space < Nh/(W 20 atom cm -2 ) < 20, < 
Z/Z @ < 2, and 1 < T/keV < 10, where N H is the column 
density of the absorbing screen, and Z and T are the metallicity and 
temperature of the cluster. Projections of the ratios for constant Z or 
Nh are shown in Fig.[Tl|. The colours were calculated for an object 
of redshift of 0.01. We found that the choice of redshift was not 
significant, due to the low redshift nature of the clusters examined. 
The predictions of the model show some differences from those of 
Allen & Fabian (1997), particularly at low temperatures, where the 
colour ratios are higher than those predicted before. 

We calculated, too, the colour ratios expected for a cooling 
flow model (Johnstone et al. 1992) with an absorbing screen (Fig. 
[L?] ). We modelled a gas cooling to 0.001 keV over the same range 
of abundance, upper temperature and absorbing column density as 
the isothermal gas model. 

Our first step was to fit the C/D colour ratios for the contours 
of a cluster with a single component isothermal model. We assumed 
galactic absorption, Nn,^x, for each cluster and a constant metal- 
licity of Z c \ ust = 0.3 Zq. By taking the average fitted temperature 
of the outer three contours we derived an isothermal temperature of 
each cluster, Ti so . Table ^ shows the fitted isothermal temperature 
for each cluster, with the galactic absorption used. The galactic col- 
umn densities listed were obtained using the FTOOLS NH program 
(Dickey and Lockman 1990), except for A478, PKS 0745-191 and 
A31 12, for which we used the values quoted in White (2000). 

Temperature profiles calculated using an isothermal model, 
with fixed absorption and metallicities, are shown for a sample of 
four clusters in Fig. |n| The errors in the temperatures were prop- 
agated directly from the errors in the observed C/D ratio only. 
Three of the clusters are known cooling flow clusters. Applying 
the isothermal model to each contour shows a temperature de- 
crease from the outer regions of the clusters to their centres. For 
the non-cooling-flow cluster there is no evidence of this tempera- 
ture change. This is generally the case for our sample of clusters. 
Those with strong cooling flows show a temperature decrease, ex- 
cept for those masked by a strong galactic column. 



4.1 Fitting metallicity gradients 

The cooling flow clusters show a temperature decrease from their 
outer regions to their centres, assuming a single-phase isothermal 
gas. White (2000) found that 90 per cent of his sample of 98 clus- 
ters were consistent with isothermality (after including cooling flow 
effects) at the 3-cr level. We attempted to fit the contour colours 
with an isothermal gas with a varying metallicity. We minimized 
the function 



x 2 (z) = 



f [B/D] ohs - [B/D] iao (Z,T iao ,N H ^i) 

I a lB/D] obs 



+ f [C/D] ODS - [C/D] iso (Z, T iso , N H , sal ) \ 2 (4) 

I °"[C/D] oba J 



to find the optimum value of Z, where [B/D] h a and [C/D] bs 
are the observer contour colour values, c"[B/D] obs and o"[c/D] oba 
are their respective errors, and [B/D]i so and [C/D]i so are the pre- 
dicted isothermal values as a function of Z. 

This analysis showed that simply using an isothermal gas with 
a metallicity gradient was not a good fit to the observed data in the 
inner parts of the clusters. The fitted metallicity was 0.9Zq for the 
inner contour of the Virgo cluster data, but the \ 2 value was ~ 600. 
Most of the other strong cooling flows showed similar results. 



4.2 Fitting cooling flow models 

To fit our data it was clear that another component to the model 
was required. We added an extra cooling flow component, which 
was an obvious first candidate. If the number of counts observed 
in a particular band is denoted by X b B , an isothermal model pre- 
dicts a count of X- mo {Z, T, Nh), a cooling flow model predicts 
X C {(Z,T, Nh), where T in this case denotes the upper cooling 
temperature, and the models are normalised by the values N[ so and 
Net, then we can write 

Bobs = N iso B iso (Z,T,N H ) + N c{ B cf (Z,T,N H ), (5) 
Cobs = N iso C iso (Z, T, N H ) + N c t C c! {Z, T, N H ), and (6) 
Dobs = N iso D iso (Z,T,N H )+N ci D cf (Z,T,N H ). (7) 

By solving the above equations we can find the cooling flow 
fraction in a particular band. In Band C the cooling flow fraction is 



fo 



N cl C cf 



Ni so Ciso + Nd Ccf 

C c f(C is o - [§] obs Aso) 



(D c f Ciso — C c f -Diso) 



(8) 



Note that the model predictions, Xi so and X c f are functions of 
metallicity, absorbing column density and temperature. We made 
the assumption of constant metallicity, isothermal temperature (or 
upper cooling temperature) and obscuration by galactic column. 
Those clusters with strong cooling flows show large cooling flow 
fractions at their centres. In Fig. M4| we present C band cooling flow 
fraction profiles for our four example clusters. The errors shown 
only take account of the uncertainties in the C/D ratio. 



4.3 Calculating intrinsic absorption 

By rearranging equations <Jsj> , (^) and (^), we can predict the ob- 
served B/D colour from the C/D colour and the models. Rear- 
ranging, 



'.%] 



Ccf 



prccl 



B iso D cl 



+ 



-1 

- D J obs 



Diso 



C i30 Dr 

B,. ( D,„ 



(9) 



We found that the predicted B / D colour does not match the 
observed value. The B/D ratio is more sensitive to absorption than 
the C/D ratio (Fig.[Tl| and[I|). Assuming the difference between 
the observed and predicted values of B/D is due to intrinsic ab- 
sorption, then it is possible to use the above equation to find the 
increase in absorption. We varied the obscuring column density to 
find {B, C, -D}{iso.cf} until the predicted B/D ratio from equation 
^ was the same as the observed one. 

We then plotted a profile showing the required extra absorp- 
tion above galactic as a function of contour for each cluster. Those 
clusters with strona cooling flows show an increase in absorption at 
their centres. Fig. |^| shows profiles of the required extra absorption 
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for our example four clusters. Listed in Tableware the required ab- 
sorption increases for each cluster in our sample. This value is the 
difference between the absorption in the central contour and the 
mean of the two outer ones. We took the mean outer contours to de- 
crease the chance of background contamination. We looked at the 
difference in absorption between the centre and outside of the clus- 
ter, rather than the absolute difference from the central absorption 
to galactic, because the galactic absorption value is uncertain for 
many clusters. The errors on the values assume that all uncertainty 
is due to the observed B/D ratio. S(ANh) shows the significance 
of the increased absorption using an expression with the same form 
as equation (0). 



5 DISCUSSION 

5.1 Adaptively-smoothed images 

Many of the cooling flow clusters in our sample show relatively 
featureless elliptical images after adaptive smoothing. We present 
some interesting ones, and those of the Perseus, Virgo and A2199 
clusters in Fig. |l| ^| and ^| The clusters show structure which ap- 
pears to vary between the different bands. 



effect (Fig. |ll|), so if we have overestimated the levels of absorp- 
tion, then it will be primarily for those clusters with metallicity 
gradients. However, we found that metallicity gradients were not 
sufficient by themselves to account for the observed colour gradi- 
ents. 

The PSPC detector limits the energy bounds of our observa- 
tions. There are too few independent energy bands to make more 
than two colours, which limits the number of physical quantities 
we can fit the data. More colours will help to completely resolve 
any degeneracy between metallicity and absorbing column density. 
Data from Chandra will be able to solve this completely. The in- 
creased effective area of the telescope will also improve our statis- 
tics, reducing the uncertainties in our colours. 

Allen & Fabian (1997) used a partial screening model in con- 
junction with a deprojection method on their data. Using a partial 
screening model results in larger column densities than those pre- 
dicted without. Doing a similar analysis would show larger levels of 
absorption in our clusters, however it would add unnecessary com- 
plexity to our procedure. We also did not take the absorption results 
to check for consistency with the C / D ratio, but do not expect the 
results to be significantly different. 



5.2 The sample of clusters 

Fig. ^ and|^ show that clusters with low galactic latitude, or those 
masked by regions of high absorption, show low B/D and C/D 
colours. These clusters include, for example, Ophiuchus, PKS 
0745-191 andCygnus A. 

Those clusters with strong cooling flows show large C / D ra- 
tios, or softer X-rays, at their centres, with large C/D gradients 
moving towards their centres. These include, for example, Virgo, 
Perseus, 2A 0335+096, Centaurus and A2199. They therefore con- 
tain cool gas at small radii, within a few arcminutes of the centre. 
This can be seen by comparison with the theoretical curves in Fig. 
[n[ The only way of achieving the gradient is by adding cooler gas. 
Those clusters without strong cooling flows do not show this strong 
colour gradient, for example, Coma, A2256 and Al 19. 



5.3 Comparison with models 

The isothermal temperatures derived using the isothermal model 
in the outer regions of the clusters agree in most cases with those 
derived from ASCA data (White 2000). There are a few clusters, 
most of which are absorbed by high galactic column densities, for 
which we find unusually low temperatures. These include Perseus, 
PKS 0745-191, A478 and Cygnus-A. We note that we were able to 
reproduce the Perseus ASCA temperature result using a column of 
10 21 atom cm" 2 . 

Those clusters with strong cooling flows show good evidence 
of absorbing material at their centres. For example Virgo, Centau- 
rus, Perseus and A478 all show significant increased levels of ab- 
sorption. There is no evidence for absorbing material in those clus- 
ters without cooling flows. Abundance gradients in cooling flow 
clusters have been observed, for example in the Centaurus cluster 
(Fukazawa et al. 1994), the Virgo cluster (Matsumoto et al. 1996) 
and AWM7 (Ezawa et al. 1997), with the metallicities decreasing 
outwards from the cluster centres. Due to the relationship on B/D 
colour from absorption and metallicity, there is some degeneracy in 
the two variables. At low temperatures metallicity has its strongest 



6 CONCLUSIONS 

We have presented an analysis of X-ray colour maps. The analysis 
technique attempts to group areas of the image together in order to 
maximize the signal to noise of the results, but preserving informa- 
tion about the structure of the object. This technique will be useful 
in examining data from the current generation of new X-ray tele- 
scopes. These data will show large variations in count rate, and will 
need analysis which is more intelligent than simple binning. 

We created colour X-ray maps of the cores of a sample of 33 
clusters, almost doubling the sample of Allen & Fabian (1997). The 
profiles generated from these maps show that there is cooling gas 
at the centres of those clusters with strong cooling flows. There is 
no evidence for cooling gas in clusters inferred from imaging de- 
projection or other methods to have weak, or no, cooling flow. We 
find a clear anti-correlation between the central cooling time of a 
cluster and the C/D colour gradient significance. The central cool- 
ing time is an indicator of the quality of cooling flow data. We also 
find a correlation between the cooling flow flux and the C/D sig- 
nificance. The maps of the non-cooling-flow clusters contain more 
structure as a group than those clusters with cooling flows. 

We fitted a single-phase isothermal model to the C/D colour 
of each of the contours of our sample of clusters. Those clusters 
with strong cooling flows show a significant decrease in the fitted 
temperatures of their inner contours relative to their outer contours. 
An isothermal gas with a varying metallicity alone was also not able 
to fit our data. We added a cooling flow component to the model 
with the same upper temperature as the isothermal plasma. Those 
clusters with cooling flows required the addition of a significant 
fraction of cooling flow model in their centres to account for the 
C / D colour. To fit the observed B/D colour we required the addi- 
tion of extra absorbing material at the centres of our strong cooling 
flow clusters. This provides more evidence that cooling flows accu- 
mulate cold material. There was no evidence for intrinsic absorbing 
material for clusters without cooling flows. 
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